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I. INTRODUCTION 



Topological defects appearing in different field-theoretical models and residing in higher- 
dimensional space-time can be considered as a prototype of four-dimensional world provided 
the known particles and gravity are localized on them ||^. For example, a domain 
wall solution of a simple theory with spontaneous breaking of ^ —cf) symmetry in 
five-dimensional space-time leads naturally to four-dimensional chiral fermions residing on a 
domain wall and to four- dimensional gravity if some fine-tuning of a bulk cosmological 
constant to the domain wall tension is made. A number of "thick" wall solutions in scalar 
field theory coupled to gravity has been found and the possibility of gravity localization 
on these walls has been studied. 

If a domain wall is replaced by a more complicated topological defect, such as string 
|[TT| , monopole [|l2|,P,p!3| or instanton |I^JT5[] (living, respectively, in 6, 7 or 8 dimensions), the 
structure of chiral fermionic zero modes, needed for construction of a realistic phenomenology 
in four dimensions, gets richer. Hence, if the possibility of constructing standard model 
interactions along these lines from higher dimensions is taken seriously, other fields (scalar, 
gauge and gravity) must be localized on topological defects as well. 



In was shown that a thin string in 6 dimensions could lead to localization of gravity if 
certain relations between the tension components are satisfied (for even higher-dimensional 



generalisations see IH^H)- The aim of the present paper is to provide a field-theoretical 
realisation of this idea and show that an Abelian Higgs model Jl^, in which a string defect 
arises naturally, can lead, at once, to solutions which localize gravity and which have a 
geometry free of singularities. In the absence of branes a model with an Abelian magnetic 
field in six dimensions was considered in [Q. 

The plan of our paper is the following. In Section II the Abelian-Higgs model coupled 
to six-dimensional Einstein-Hilbert gravity is formulated. In Section III the properties of 
six-dimensional cylindrically symmetric warped geometries will be discussed in the absence 
of defects. This is a necessary step for a classification of the solutions in the presence of 
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vortex-type defects. In Section IV the asymptotics of the solutions will be studied both 
for the metric and for the scalar and gauge fields. Necessary conditions in order to get 
string solutions will be formulated. Section VI contains numerical examples of the solutions 
localizing gravity whose parameter space is analyzed in detail. The thin string limit is 
considered. In Section VII we discuss solutions leading to gravity localization in the case of 
higher winding numbers and consider regular string solutions which do not localize gravity. 
Section VIII contains our concluding remarks. In the Appendix useful results concerning 
the interplay between the thin string limit and the Bogomolnyi limit have been collected. 



II. ABELIAN-HIGGS MODEL IN SIX DIMENSIONS 
A. Basic equations 

The total action of gravitating Abelian Higgs model in six dimensions can be written as 

where 5'brane is the gauge-Higgs action and S'grav is the six-dimensional generalization of 
Einstein-Hilbert gravity. More specifically, 

^brane = j ci'xv^^brane, C^r.ne = \{V ACf))* cf) - ^FabF^'' - ^ (0*0 ' V^)\ (2.2) 

where Va = Va — ^^^a is the gauge covariant derivative, while is the generally covariant 
derivative [|. In Eq. ( p. 21) , v is the vacuum expectation value of the Higgs field determining 
the masses of the Higgs and of the gauge boson 



rriH = y2Xv, my = ev. (2.3) 



The conventions of the present paper are the following : the signature of the metric is mostly 
minus, Latin (uppercase) indices run over the (4 + 2)-dimensional space whereas Greek indices run 
over the four-dimensional space. 
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As usual the bulk action is 

~ R 



Sgra.v — ~ I d^Xy—G 



+ A 
2X 



(2.4) 



where A is a bulk cosmological constant, x = SvrGe = Sh/Mq and Mg denotes the six- 
dimensional Planck mass. 

From Eqs. (|2.1| )- (|2.4| ) the coupled system of classical equations of motion can be derived 

G^-^VaVb0 - e^AAA'^(j) - ieAAd^cj) - leW A^A^ct)) + A(0*0 - = , (2.5) 

y^F^B^_^2^B^*^^^_^(^^QB^*_^*QB^'^ ^ ^2.6) 

Rab - \GabR = X (Tab + AGab) , (2.7) 

where 

Tab = -jC^r.ncGAB + I [iVA<pyVB<p + iVB<pyVA<p] - FacFb^ . (2.8) 

We are interested in static string-like solutions of Eqs. ( p.5| )-( pr7D that depend only 
on the extra coordinates and thus do not break general covariance along the four physical 



dimensions. For these solutions a six- dimensional metric is POl,!!! 



ds^ = GABdx^dx^ = M^{p)g^^dx^'dx'' - dp^ - L{pfd9\ (2.9) 

where p and 9 are, respectively, the bulk radius and the bulk angle, g^^ is the four- 
dimensional metric and M{p), L{p) are the warp factors. The Nielsen-Olesen ansatz for 
the gauge-Higgs system reads: 

<jy{p,e)=vf{p)e'^', 

Aeip,9) = -Jn- Pip)], (2.10) 

where n is the winding number. 

Inserting Eq. ( |2.1C ) into Eqs. ( |2.5| )-(2!7) and using the specific form of the metric given 



in Eq. (p^.9|), the equations of motion can be written as 
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/" + (4m + i)f + (i-f)f-^f = 0, (2.11) 

P" + {Am- i)P' - afP = 0, (2.12) 

i' + 3m' + f + 6m^ + 3im = -/i - z/tq + (2. 13) 

4m' + lOm^ = -/i - z/re + (2.14) 

2m£ + 3m^ = -|-^r, + ^. (2.15) 

The four-dimensional metric g^u{,x,t) obeys the equation: 

- = ^Aph^,. , (2.16) 

where the physical cosmological term, Aph, is an arbitrary integration constant and Mp is 
the four-dimensional Planck mass. If Aph = we can choose g^i, = rj^y where rj^^ is the 
Minkowski metric. This will be the case we will be mostly interested in, even though we 
will often keep the physical cosmological constant, for sake of completeness. 

In Eqs. ( p.ll|) -( p.l5|) the prime denotes the derivation with respect to the rescaled 
variable 

x='^. (2.17) 
The function L{p) appearing in the line element of Eq. (|2.9| ) has also been rescaled, namely 



^(^) = ^^(P)- (2.18) 

Furthermore, the functions m(x) and [{x) are simply 

d\TiM{x) d\nC{x) 

= ; , t{X) = ; . (2.19) 

ax ax 

In Eqs. (|2.11| )- (|2.15| ) dimensionless parameters have been defined as: 



rriy 2xA mjr IGvrAph 

a = 2—, /i = u = XiTT-, /iph = , .2 2 • (2-20) 



Within the mentioned conventions the components of the energy-momentum tensor 
turn out to be 



r„(.) . T» = i;' = + \(f - 1)^ H- ^ + ^, (2.21) 
r,W.T; = -^ + -(f- 1)^-^ + 1^. (2.22) 

/•/2 1 p/2 f2p2 

r„(x).T;=^ + j(f- 1)^-^-1^. (2.23) 

Eqs. ( |2.11|) -( ^.12|) correspond, respectively, to the equations of motion for the scalar and 
gauge fields whereas Eqs. (|2.13| )-( ]2.15|) correspond to the various components of Einstein's 
equations. Eqs. ( p.l3|) -( ^.14| ) are truly dynamical since they involve m!{x) and l'{x). Eq. 
( p.l5| ) connects m{x) and ^{x) to the first derivatives of the gauge and scalar degrees of 
freedom and it is, therefore, a constraint. In the following part of the paper /Xph = will be 
assumed. 



B. Boundary conditions 

Regular solutions of the equations of motion can be investigated once the boundary 
conditions are specified. In order to describe a string-like defect in six dimensions we have 
to demand that the scalar field reaches, for large p, its vacuum expectation value, namely 
V for p oo. In the same limit, the magnetic field should go to zero. Moreover, 
close to the core of the string both fields should be regular. These requirements can be 
translated in terms of our rescaled variables as 

/(O) = 0, hm f{x) = 1, 

P(0) = n, lim P{x) = 0. (2.24) 

The requirement of regular geometry in the core of the string reads: 

M'(0) = , £(0) = , £'(0) = 1 , (2.25) 

and one can arbitrarily fix M(0) = 1. 

These conditions listed in Eqs. ( |2.24| )-( |2.25| ) completely specify a solution. Clearly, Eqs. 



( p.24| )- (|2.25| ) do not tell anything about the asymptotics of the metric at large distances 



from the string and thus a possibihty of gravity locahzation on a defect remains open. The 
requirement of gravity locahzation is equivalent to the finiteness of the four dimensional 
Planck mass, 

Ml = f dxM^(x)C(x) < oo , (2.26) 

mjj J 

which does not hold in general and requires a fine-tuning of parameters of the model, as it 
will be discussed. 



III. BULK SOLUTIONS 



Outside the core of the string all source terms in Eqs. ( 2.21|) -( ^.23 ) vanish and a general 
solution to equations 



dm 5 n a 
/i 3m 



4m 2 ' ^^'^"^ 



can be easily found 20 1. We will consider the case A < only, since the case of a positive 



cosmological constant in the bulk was studied in Defining 



> 0, (3.2) 
10 ' ^ ^ 



the solutions of Eqs. ( p.l| ) can be written as 



where e is an integration constant. Once m{x) is known, £{x) can be determined from the 
second of Eqs. (pTip. 

In order to understand whether gravity can be localized on any of those solutions we 
compute M and C: 



M{x) = Moe-'^\l + ee" 

1 1 + ee^^^' 1 1 



Cix)=Coe~'''^ (3.4) 



If e > or e < — 1 we have that for x — oo M ~ £ ~ e*^^. For this solution the integral 
in ( |2.26| ) diverges and thus gravity cannot be localized. 



If e = the solution is simply 

m(x) = i{x) = — c, (3.5) 
and the warp factors will be exponentially decreasing as a function of the bulk radius: 

M{x) = Moe-"^ 

jC{x) = ^oe"'^ (3.6) 

The solution of Eqs. ( p.5| )-( p^ leads to gravity localization and to a smooth AdS geometry 
far from the string core. We also note that a pure positive exponential solution can be 
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derived for e ±oo taking Mq ~ Ci/|e|5 and Cq ~ C2/|e|5 where Ci and C2 are two arbitrary 
constants. 

If — 1 < e < 0, M(xo) = where Xq = logl/|e|. The geometry is singular at x = Xq 
(see below) and x < Xq should be required. In spite of the fact that C diverges ai x = Xq, 
the integral (|2.26|) defining the four dimensional Planck mass is finite. So, these geometries 
can be potentially used for a gravity localization, provided the singularity at Xq is resolved 
by some way, e.g. by string theory [|19|. These types of singular solutions were discussed in 
pO| for the case of positive cosmological constant. 

Finally, if the bulk cosmological constant is zero, the solutions have a power-law be- 
haviour, namely 

M{x) ~ x^, C{x) ~ x\ (3.7) 

with 

^7 + (5 = 1, d-f^ + S^ = 1, (3.8) 

where d = 4 is the number of dimensions of the metric Qf^^. These solutions belong to the 
Kasner class. The Kasner conditions leave open only two possibility: either 6 = 1 and 7 = 
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or 7 = 2/5 and S = —3/5, [the same fractions as in Eqs. ( p.4| )]. None of them could lead to 
localization of gravity. 

In order to study the singularity structure of the bulk solutions without specifying e, the 
different curvature invariants may be computed. In the case of warped metrics the Riemann 
invariant has been calculated in |^. The explicit form of all the curvature invariants for 



the metric of Eq. (|2.9|) is given in [|l9l and they are: 



R = 20m2 + 8m + 2i' + 21^ + Sim, (3.9) 
RabR^^ = 80m^ + 20m'2 + 2i'^ + 2t + Unirn^ + Al'^l' + 28m^^^ 

+ 32m^£ + 8mm7 + 8m7' + Sm'^^ + Sm^f + SmW + 8m£^ (3. 10) 

RabcdR^^^^ = - 24m^ + Sfi' + 4f ^ + 32m'^m' + 16m'^ (3.11) 

CAiJCD^^^^^ = ^ [{m' - t) + i{m - . (3.12) 



Let us first check the regularity properties of the bulk solutions (|3.3| ). The Ricci and 
scalar curvature invariants are simply constant for any e: 

R = 30c2, RabR^^ = 150c^ (3.13) 

The Riemann and Weyl invariants are 

nABCD _ -20 (e-2°"^ - 12 e-i^^^ e - 138 e"i°=^ - 12 e'^"^ + e^) 
RabcdR , ,4 , 

{pi — OCX I - ' 



3840c^e-^°^^e2 



CabcdC^'''''' = ^-^^f—^ (3.14) 

From these last two expressions we can see that ife>Oore<— 1 the Weyl and Riemann 
invariants do not have any pole and thus we have a regular geometry. For e = all curvature 
invariants are constant and we have a six-dimensional AdS space. On the other hand if 
— 1 < e < both Weyl and Riemann invariants diverge at a finite value of x leading to a 
singular geometry. The case e = — 1 is somewhat specific since a singularity of Kasner type 
is developed in the origin. 

The parameters e, Mq, £o are not fixed for a bulk space solution. If a string-like defects 
is placed at the origin x = these constants are no longer arbitrary and are functions of 
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the three parameters of the model, namely, e = e{a, fi, z/), and so on. An everywhere regular 
geometry is simultaneously achieved together with gravity localization if the parameters lie 
on the surface e{a, /i, z/) = 0. 



IV. ASYMPTOTICS OF THE SOLUTIONS 

A. The asymptotics of the solutions at the origin 

The form of the solutions in the vicinity of the core of the vortex can be studied by 
expressing the metric functions together with the scalar and gauge fields as a power series 
in X, the dimensionless bulk radius. The power series will then be inserted into Eqs. (p.ll )- 



( P.15| ). Requiring that the series obeys, for x ^ 0, the boundary conditions of Eqs. (p.24|) 



the form of the solutions can be determined as a function of the parameters of the model. 

Consider, in particular, the case of = 1 0. In this case asymptotic solutions close to 
the core can be written as : 

+ - + — — + — 1 + 25 -x^ (4.1) 



3 6 3a 3 
P{x) c^l+Bx^, (4.2) 



C(x) ~ a; + 



/i , 1 



x-\ (4.4) 



12 M8 Qa 6 

In Eq. ( [4.4|) A and B are two arbitrary constants which cannot be determined by the local 
analysis of the equations of motion. These constants are to be found by boundary conditions 
for f{x) and P{x) at infinity. 

For completeness, we give also the limit of the components of the energy-momentum 
tensor in the vicinity of the origin: 



^This analysis can be generalized to the case of higher winding number (namely n > 2) as it will 
be seen in Section VI. 
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To 



I p — 



re 



- + + 

4 a 
1 2fi2 

4 a J 
1 



4 



(4.5) 



B. The asymptotics of the solutions at infinity 



Consider now asymptotics of solutions at large x, assuming that the geometry is regular 



at infinity. This situation can be realized both in the case e = [i.e. Eq. (|3.5| )] and in the 
case Eq. ( p.4|) either with e > or with e < — 1. We can write 



(4.6) 



P[x) = P + 5P{x), 
f{x)=J-6f{x), 

where, according to Eqs. ( p.24| ) for x 1, P ~ and / ~ 1. 

First, we will consider the case e = 0. From Eq. ( |2.12D it follows that 

5P" - 3c6P' - a6P = 0, 

which implies that, for a; ^ 1 



(4.7) 



(4.8) 



If 4a ^ 9c^ (the limit of small bulk cosmological constant) the solution is compatible with 
the gauge field decreasing asymptotically as SP ~ e~^^ which is the well known fiat-space 
result of the Abelian-Higgs model. 

The linearization of the scalar field equation, using a similar procedure, yields 



Sf - 5c6r - 26 f = 0, 



(4.9) 



leading to the solution 



5c, 



(4.10) 
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If 25c^ <^ 8 the perturbed solution goes as 5/ ~ e which, again, is the flat-space result. 
In Eqs. ( |4.8| )-( ^?TI1 ) a minus sign should be chosen to obtain an exponentially decreasing 



behaviour. Moreover, linear approximation is only valid if o"2 > 2cri + 2c, in the opposite 
case bf ~ exp(2cri + 2c). 

For e 7^ the sign in front of c in equations (|4.7| , [4.9|) must be changed, and we get 



The linear approximation is valid if 02 > 2ai — 2c, otherwise 6f ~ exp(2cri — 2c). 



C. String tensions 

In the description of string-like defects in four-dimensions it is useful to study the relations 
obeyed by the string tensions. Indeed, some particular linear combinations of the string 
tensions may be related to the Tolman mass per unit length and to the angular deficit of 
the geometry. These quantities can be used in order to classify the physical properties of 
the solutions and in order to check if they are effectively asymptotically conical P2 . 



This analysis is also relevant in our case even though the physical interpretation changes 
and the specific algebraic relations are different from the well known four- dimensional cases. 
The components of the string tension are defined as 

poo 

fii= dxM^{x)C{x)Ti{x). (4.12) 
Jq 

The integrals appearing at the right hand side always converge for solutions with regular 
geometry, as follows from the results of the previous subsection. 

We will show now that for M(x) — > 0, £i{x) at a; ^ 00 (i.e. in the case we are 
mostly interested in) the constant B in Eq. ( [4.4[ ) can be analytically determined. 

Consider a specific linear combination of Eqs. ( |2.13| )-( p.l5D , namely 
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2m! + 4m^ + ml = - 
e' + f + 4im = 



V 



1 3 
^0 + -Jp - -Te 



4 



4 



(4.13) 
(4.14) 



Integrating Eqs. ( |4.13| )-( |4.14| ) from zero to Xc — >■ oo, we get 



M{x,yC\x 



1 - 



-J M Cdx - l^/io - -/ie + — 



(4.15) 
(4.16) 



In the limit Xc ^ oo, Eq. ( [4.15|) is the six-dimensional analog of the relation determining 
the Tolman mass whereas Eq. ([4.16|) is the generalization of the relation giving the angular 
deficit. Therefore the following relation must hold: 



fJ'O ~ l^e — — 

V 



(4.17) 



This is the fine-tuning relation which has been already found in [|TT],|T2|. 

Let us note that if e 7^ one can use the asymptotics of the metric found in Eq.(p.4|) for 
computation of the left-hand sides of ( [4.15| , ^?TBD . This gives a more general relation, 



/io - /ie = -(1 - lOceMg £0) • 



(4.18) 



The extra term vanishes also for e ±00 given the scaling behaviour of Mg ~ |e| 5 and 



~ kl^^- So, the condition ( [4.17| ) is only necessary but not sufficient in order to have 
solutions leading to gravity localization. 

Being an integral relation we can use ( |4.17] ) in order to get informations on the behaviour 
of the solution in the origin. In order to do so, the left hand side of Eq. ( |4.17| ) can be directly 
expressed in terms of the scalar and vector fields 



p/2 f'^P'^ 

+ 



dx 



(4.19) 



Using now Eq. ( 2.12| ) we arrive at 



C 



a- 



pPM^ 



(4.20) 
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which can be inserted back into Eq. ( [4.19|) . Integrating by parts the obtained relation the 
term containing P'^ cancels and we arrive at 



IJ'O- IJ-e = - 
a 



1 rPP' 



P P' 



(4.21) 



For the solutions we are interested in the boundary term at infinity vanishes. Using now 



the asymptotic behaviour for small x of Eqs. ( [4.4|) we obtain, from Eq. ([4.21 



1 P' 
a C 



(4.22) 



In this case, from Eq. ( |4.17|) and using Eq. ([4.21| ) we have that 

V P' _ ^ 
a C 



According to Eq. (O), for x ^ 0, P ~ 1 + Bx^. Using Eq. 
be exactly computed 



B = - 



a 
2^ 



(4.23) 

the expression for B can 
(4.24) 



This relation is very helpful for numerical integration of equations of motion performed in 
the following section. 

Finally, we would like to notice that the magnetic field can be expressed as B{p) = 



-P'{p)/{eL{p)). Using now Eq. {]^^ we have that 87r;B(0) = eM^ 



In closing the present Section we want to stress that the results obtained from the 
analysis of the relations among the tensions are only necessary conditions in order to find 
e = solutions. They are not sufficient because /'(O) ~ A still remains undetermined. 



V. SOLUTIONS WITH LOCALIZATION OF GRAVITY 

In this Section, a numerical search for solutions that localize gravity will be performed. 
The general method of integration will be firstly outlined. Secondly, some interesting nu- 
merical solutions will be presented as an example. Finally, the scan of the parameter space 
of the model leading to e(a, /i, z/) = will be presented. 
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A. Numerical integration 



In order to prepare the ground for the numerical integration the system of Eqs. ( |2.11| )- 
( ^.15| ) can be expressed in a form of a first order non-hnear dynamical system. By linearly 
combining Eqs. (2.1l)-( P^.15D the following set of equations can be obtained: 





9, 




P' = 


Q, 




C' = 


I, 




9 = 


- ^4m + 




Q' = 


— ^4m — 




m! = 


5 2 

— m — 
2 


4 



p2 
£2' 



(5.1) 
(5.2) 



M' = mM. (5.5) 



In order to solve this system the shooting method |^ has been employed. In the shooting 
method a boundary value problem is tackled as a series of Cauchy problem. Using forward 
integration in x with a specific set of initial conditions, the asymptotic values of the fields are 



found. If the asymptotic values are different from the boundaries discussed in Eq. ( p. 241) , 



the integration is performed again with a different set of initial conditions until the correct 
asymptotic boundary is reached. 

The technical problem is twofold. On one hand the shooting method requires that the 
initial conditions are given in a range of the parameter space sufficiently close to the one 
admitting e(Q;,yU, z/) = solutions. On the other hand, for the scalar and gauge fields the 
integration is a boundary value problem whereas, for the warp factors, the integration is 
a Cauchy problem. The shooting method is implemented in such a way that each forward 
integration is performed by a Runge-Kutta routine. The differences between the obtained 
values at large x and the boundary values of Eq. ( p. 241) are considered as functions of the 
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constants A and B appearing in the asymptotic solution in the vicinity of the core. The 



Newton method I^S] is then used in order to find the values of A and B required in order 
to match the boundaries within the precision of the algorithm. 

Therefore, in order to identify a correct set of initial conditions extensive use of back- 
ward integrations has been made. Backward integrations have been applied by giving final 
conditions on warped solutions at large x and by integrating the equations from large to 
small X. This method allowed to identify a set of initial conditions leading to the correct 
asymptotic boundaries with warped metrics. This idea has also been applied in a different 
context always involving warped compactification in string inspired models with quadratic 
curvature corrections [1191. 



B. Some numerical examples 

We are now going to give some examples of the numerical integration, see Fig. |T]. The 
values of the scalar and gauge fields are the ones dictated by Eqs. ( p.24| ). At the core 
of the vortex, the behaviour of the metric and of the other fields follows Eqs. ( [4.4| ). For 
large x the metric functions decay exponentially. The asymptotic behaviour of the metric 
functions for large x belongs to the class of solutions discussed in Section II for e = 0, namely 

m i — a/ — /i/10. 

In the example of Fig. |l|, a = 1.16. The relationship ([4.24|) is satisfied. The curvature 
invariants (scalar curvature, Riemann, Weyl and Ricci invariants) defined in Eqs. (|3.9|) - 
( p.l2| ) are plotted in Fig. ^ for the solution illustrated in Fig. |l|. The Riemann and Ricci 
invariants go asymptotically to a constant which is 0(c'^), whereas the scalar curvature goes 
to a constant 0{c^) as dictated by the asymptotic solutions with e = and derived in Eq. 
( p. 51 ). The Weyl invariant goes to zero since, asymptotically, the space-time is isotropic. 
Finally it is interesting to check the behaviour of the various components of the energy 
momentum tensor, which are reported in Fig. ^. From Fig. ^ we notice that Te{x) gets 
negative and it approaches zero from negative values. This is due to the fact that, in Tg{x) 
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a=1. 1600000000 v=2.5801 288700 n=-0. 0782271 150 
A=0.481 1726769749937183 B=-0.2247949774303944492 




FIG. 1. A vortex-like solution of the Abelian-Higgs model leading to gravity localization. With 
the full thick lines we report the behaviour of the scalar and gauge fields whereas with the thin 
lines the behaviour of the warp factors is illustrated. The same convention will be used throughout 
the paper. 



R 

R Rl"^ 

fivctfS 
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15 



FIG. 2. The curvature invariants for the solution shown in Fig. ^. 
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FIG. 3. The components of the energy momentum tensor of the sources in the case of the 
numerical solution described in Fig. ||. 

the magnetic stresses dominate leading to locally negative pressure. 

In Fig. H, we show an example where the exponential damping of the metric is slower 
than in the case showed in Fig. |l]. The reason for this behaviour stems from the fact that 
the value of the cosmological constant (i.e. the /i parameter) chosen in Fig. ^ is roughly one 
order of magnitude smaller than in the case discussed in Fig. ^ The case a = 2 corresponds 
to mv = rriH, i.e. Bogomolnyi limit where the Abelian-Higgs model in four dimensional flat 
space obeys interesting properties. This limit is discussed in some details in the Appendix. 
To check for the accuracy of numerical integration procedure we consider the constraint 
equation (|2.15|) (which was not used in the numerical procedure) that connects m(x) and 
£{x) to the Tp{x) component of the energy-momentum tensor: 

2mi + 3m^ + ^ + -Tp = 0. (5.6) 

Inserting the numerical solutions obtained in the present Section into Eq. ( ^.6|) we should 
find that Eq. ( ^.6|) is numerically satisfied. The accuracy for which Eq. (|5.6| ) is satisfied 
indicates the accuracy of the integration. 
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a=2.0 v=2.157 |i= -0.0035408066 
A= 0.582760000786 B=-0.46360369837371 
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FIG. 4. A warped solution is illustrated in the case a = 2 where mn = my- 

In Fig. ^ we plot the left hand side of Eq. (|5.6| ) for the solution shown in Fig. |l] of this 
Section. As we can see Eq. (|5.6| ) is satisfied with a precision of 10~^^ in the vicinity of the 
core and with a precision of 10^^^ for large x. For the other solutions similar behaviour can 
be obtained for the constraint. 

C. The parameter space of the solutions with gravity localization 

Up to this moment, examples of compactifications on a vortex were analyzed. A scan in 
the parameter space of the solution will now be performed. 

In order to achieve this program in a systematic way the procedure is, in short, the 
following. Suppose that we have already a solution that satisfies all the requirements (for 
example, the one in Fig. |l] or It can be found by a procedure of backward integration 
discussed above. Then one of the three parameters among a, v and /i is kept fixed, a second 
is changed by small steps and the third one is tuned to the first two by requiring that 
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a=1 .1600000000 v=2.5801 288700 |i=-0.0782271 150 
A=0.481 1 7267697499371 83 B=-0.2247949774303944492 




FIG. 5. The constraint Eq. (5^) is illustrated for the case of the solution given in Fig. 



{x) + 



10 



10 



(5.7) 



are satisfied for large x. The zeroes of Eq. ( p. 7] ) are found by means of a simple bisection 



method 23 



Using this procedure a scan of the parameter space of the model has been performed. 
Let us start by discussing the scan in the (/z, a) plane at constant u. The results of the 
study are shown in Fig. ||. 

The same type of procedure can be carried on in the case of fixed fi. As a result we will 
get the (a, z/) projection of the parameter space of the solution. This result is presented in 
Fig. 1^. Finally we can also investigate the (/x, z/) plane at for fixed a. In this case the results 
of our study are shown in Fig. In Fig. ^ we can notice that the limit yU — > at a = 2 
implies that u 2. This can be shown analytically, see Appendix. 

Finally, in order to give a visual feeling of the fine-tuning surface we can present a three- 
dimensional plot in the z/) plane. This surface is illustrated in Fig. p. The points on 
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-0.25 -0.2 -0.15 -0.1 -0.05 

FIG. 6. We illustrate the variation of {a, /j,) for fixed v. Along each line in this plane the 
numerical u is constant and it is given by the numerical values reported at the bottom of the 
curves. 




1.2' 1 1 ' 1 ' ' 1 ' 1 1 

2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3 3.4 

V 

FIG. 7. The (a, v) projection of the parameter space of the warped solution is illustrated. The 
curves correspond to lines of constant /x. 
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FIG. 8. We report the variation of (/i, zv) for fixed a. As in the previous figures the numerical 
values represent the values of a which is kept constant along each curve. 

this surface defines a fine-tuned solution of the types presented in Fig. |I] and Fig. ^. 

D. The thin-string limit 



The thin string hmit was used for gravity locaUzation in |[TT|]. In this limit the energy 
density of the vortex gets more and more localized in the origin while the exponential 
damping of the warp factors gets milder as in Fig. |^. Here we will show how the thin-string 
limit implemented in the Abelian Higgs model. 

The main idea of a thin string limit is to replace the real thick string made of the scalar 
and gauge fields by a thin gravitating object Q characterized by three different string tensions 
defined in Eq. ( |4.12| ). Then the parameters of the metric outside the string Mq, Cq, e 
and four-dimensional Planck constant are related in a simple way to the string tensions 



^ The physical size of this gravitating object, Iq, shall be much smaller than the "size" of extra 
dimensions given by rn = 



10M| 
SttA • 
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-0.25 1.5 

FIG. 9. We illustrate the three-dimensional picture of the fine-tuning surface in the (a, z/) 
plane. The points on the surface correspond to warped solutions obeying the boundaries of Eq. 



(2.24). 



[i.e. /io, /ip and /ie], to the six- dimensional Planck constant, and to the bulk cosmological 



constant. These relationships were worked out for the gravity localizing solution in Ref. ||Tl 



see also they can be reproduced as well from Eqs. ([4.15|) -( fi.l6|) assuming that 



c <^ 1, which is a condition for a thin string limit in dimensionless variables. 

Here we will write somewhat more general equations valid also for e>Oore<— 1 that 
lead to regular geometry at infinity. For c ^ 1 the integrals and limits Xc ^ oo can be 
explicitly taken, and one gets 

M^Cocil - + 0{^) = ~ (/i, + lie) , (5.8) 
M^Coc{l + 8e + e^) + 0{^) = I - u (^^, - ^-^.^ + . (5.9) 

As the numerical analysis shows (see also analytic discussion in Appendix), in the limit 



c — >■ the parameters Mq and Cq go to some non-zero constants (see, e.g. Fig. |I0|). Then 



it is clear from Eqs. (|5.8|)-([5l^), that the combination /ip + fig scales as 0{c). In fact, as 
shown in [^, this feature has a more general character and it is not related to the specific 
form of the Lagrangian of the Abelian-Higgs model. 
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FIG. 10. The variation of Mq and jCq as a function of c for a = 2. 

It is interesting to ask the question whether a thin string hmit c ^ 1 can be organised 
in a way that : 

(i) the four-dimensional Planck scale is a finite constant; 

(ii) corrections to the four-dimensional Newton law stay small; 

(iii) classical gravity is applicable in the bulk; 
and, eventually 

(iv) classical gravity is applicable in the core of the string. 

In order to answer these questions let us go back to the dimension-full parameters of the 
model. They are 



A = Ml, Me, niH, A, e^. 



(5.10) 



The other parameters can be expressed as a function of the previous ones since we know 
that, for instance. 



A 



(5.11) 



H 
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In order to have only two scales in the problem, we can then take a to be fixed to given 
value (for instance a = 2). Once a has been fixed, we can see if, in the limit c <^ 1, 
the requirements (i)-(iv) can be satisfied. In order to solve the problem, we should know 
the values of Mq and Lq defined in Eq. ( p.5| ) in the thin string limit and also the limit 
of parameter v. An analytic discussion of these questions is contained in Appendix. Here 
we simply state the results. In the limit c ^ and for a = 2, z/ = 2 [see Fig. ||, and 
Mo = £o = 1 [see Fig. . 

fi) The four- dimensional Planck mass 



Ml = / dxM\x)Cix) 



m 



H 



(5.12) 



for small c and a = 2 is simply 



Mi , 

3 mH^\A.\ 

which must be equal to [1.22 x W^GeYf. 



(5.13) 



(ii) Corrections to 4d gravity were computed in |TT| and are given by ^(ro/r) . They 



are known to be small at the distances smaller than 0.2 mm |2^ and thus we require 

0.2 mm 



(5.14) 



(iii) To parameterize the Planck constant in the bulk, we write 



Mg = 7] lO^GeV . 



(5.15) 



The choice of ~ 1 makes the fundamental Planck scale to be of the order of the weak scale 



and may potentially lead to the solution of the gauge hierarchy problem [^. To satisfy 
the third requirement, we must have the curvature in the bulk to be smaller than M|, i.e. 
Ma < Mg. 

(iv) If we require the validity of gravity on the string we should demand that the 
(quadratic) curvature invariants of the solution are always much smaller than M|. In the 
small X limit the curvature invariants are 
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5(-8a + z/ {4:11 + iy)f 



64z/2 

2 



16i/2 

2 



C.BCoC--~miit^i^±^^^l^i±i^. (6.16) 



So, we find that this requirement is satisfied if 

ruH 



Me 

In numbers, the thin string limit c <^ 1 implies 



< 1 . (5.17) 



the third condition requires 

V 



c ~ 70 X < 1 , (5.18) 



I > 6 X 10-^^ (5.19) 

and the ratio ^ essential for the requirement (iv) is 

^ ^ 2 X 10-^^^ . (5.20) 

One can see that all these conditions can be satisfied in a wide range of the parameters of 
the model. For example, the hierarchy Mg ^ Mp can be easily achieved for ~ 10, ^ ~ 1. 

If the limit c ^ is taken in mathematical sense, it is impossible to satisfy all four 
conditions, because the curvature of the space inside the string diverges. However, the first 
three conditions can be satisfied, if the parameters of the model scale as: 



./1\ 4/3-3 
/I \ 2/3-1 



-2'1 



X^E-'i^^), (5.21) 
with 1 < /? < I and E being an arbitrary energy scale. 
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VI. EXTENSIONS 



A. Higher windings 

The parameter space of the solutions discussed in the previous sections has been derived 
in the case n = 1. It is also possible to discuss the case of higher windings n > 2. The 
study of these cases can be performed using the same techniques discussed in the case of 
n = 1. There are obvious differences in the asymptotics of the solutions. For the gauge field 
we have to require that -P(O) = n. This imposes a different asymptotic solution in the core, 
namely 

fix) = Ax" + (6.1) 
P(x) =n + Bx^, (6.2) 
M(x) = 1 + Cx^ (6.3) 
C{x) = x + Dx^. (6.4) 

For n > 2, we find: 

d2,, 

(6.5) 
(6.6) 

(n~5) + — — . (6.7) 



c = 


8 


V 

' 32 


+ —— 

4a 


D = 


^ + 
12 


V 

48 


6a ' 


F = 




1 + 


IBn - 



3a ' ' 6 V" 4 ^ 

With these solutions it is possible to integrate the system by repeating the same procedure 
discussed in the case n = 1. 

Since the winding number determines the value of the gauge field at the origin the 
relations among the string tensions will also be modified. In particular, Eq. ( [4.22| ) becomes, 
in the case of n > 2, 



n P' 



a C 

which implies that B, appearing in Eqs. ( |6.1|) -( |0|) , is now 

27 



(6.8) 





a=1.75 v=1 .7277073007 |i=-0.1 
A=0.1 92871 850355500001 9 B=-0.25322576331 19999740 




-0.5 L 
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15 



FIG. 11. A solution leading to gravity localization in the case of winding n = 2 is illustrated. 



B 



a 



2niJ 



(6.9) 



As in the case n = 1 the tuning of the string tensions expressed in terms of B is only necessary 
(but not sufficient) in order to get solutions which lead to gravity localization. Thus, as in 



the case n = 1, A should be again tuned. In Fig. |TT| we illustrate an example of regular 
solutions with n = 2. The components of the energy-momentum tensor corresponding to 



the solution of Fig. 11 are reported in Fig. 12 



In Fig. a solution which localizes gravity is reported in the case n = 3. In Fig. ^ 
the components of the energy-momentum tensor are reported for the solution illustrated in 



Fig. |T^. The solutions illustrated in Fig. |TT| and |T3| lead to regular geometry since all the 
curvature invariants are finite for any x. As in the case n = 1 a full scan of the parameter 
space could be performed following the same procedure outlined in the body of the paper. It 
is interesting to speculate that these solutions may have some implications in models trying 
to explain the origin of families in a higher dimensional context [^. It would then be of 
some interest to study, in the context of solutions with higher windings, the structure and 



28 



0.45 r 



0.35 - 



0.25 




0.15 



0.05 - 



-0.05 L 



10 



15 



FIG. 12. The components of the energy-momentum tensor computed from the solution of Fig 
Til are illustrated. 



a=1.75 v=1 .3080454083 (i=-0.1 
A=0. 05657531 958764981 66 B=-0.22297900731 29999939 




FIG. 13. An example of solution with n = 3. 
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FIG. 14. The components of the energy-momentum tensor in the case of the solution reported 



in Fig. 13. 



localization of the fermionic zero modes. 



B. String solutions without localization of gravity 

The examples discussed up to this point illustrated mainly solutions whose warp factors 
decrease at infinity. In this section we will give some examples of solutions with the opposite 
type of behaviour. 

In Fig. |15| we show a regular solution of Eqs. (|2.11|) - (|2.15|) where the warp factors 



increase at large distance from the core. The parameters of Fig. are a = 1, u = 1, 
H = —0.1. In the solution of Fig. |1^ Eq. (|4.24|) does not hold. From Eq. (|4.24|) we should 



have B = —0.5 whereas, for Fig. IT5| i? 7^ —0.5. In this sense the solution is not fine-tuned. 



In Fig. |1^ we plot the curvature invariants pertaining to the solution of Fig. 0. For large 
X the asymptotic behaviour illustrated in Fig. ^ can be related to the solutions given in 
Eq. 



It is interesting that the same values of a, v and /i selected for the solution of Fig. ff5 
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a=1, |j,=-0.1, v=1 
A=0.7246891 522901 01 11 29 B=-0.672927895480455777 




FIG. 17. The parameters of the solution discussed in this plot are the same of Fig. 15 but with 



different A and B. The curvature in the origin is larger than in the case of Fig. 16. 



admit a physically distinct extra solution. It is shown in Fig. and the related curvature 
invariants are illustrated in Fig. 0. We will not attempt here a full classification of the 
solutions. We can notice, however, that the solutions illustrated in Fig. ^ and Fig. can 
be distinguished (in a coordinate independent fashion) from the behaviour of the curvature 



invariants in the origin. As we can see from Figs. ^ and |lg the solution of Fig. [T^ leads to 
larger curvature invariants in the origin. 

Finally we want to demonstrate that the tuning of the string tensions as implied by Eq. 
( [4.24| ) is not a sufficient condition in order to obtain solutions with localization of gravity. In 



Fig. Oand Fig. 20 we illustrated a regular solution whose parameters are given by a = 1.75, 



u = 1.29 and /i = —0.1 and tensions fine tuned according to Eq. ([4.24|) . Nevertheless, the 



metric increases at infinity. Asymptotically for large x the solution of Fig. |T9| are the ones 
of Eq. (|3.3|) with e — > oo. 
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FIG. 18. The curvature invariants for the solution of Fig. 17 
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FIG. 19. We illustrate a solution where the string tensions have been fine-tuned according 



Eq. (4.24) but with the metric growing at infinity. 
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FIG. 20. The curvature invariants associated with the solution of Fig. 
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VII. CONCLUDING REMARKS 



In this paper we discussed the possibihty of obtaining solutions which locahze gravity in 
six-dimensions. The source of the gravitational field has been assumed to be a generalized 
version of the Nielsen-Olesen vortex living in the two-dimensional transverse space. 

We show that the localization of gravity is possible on a "thick" string and found a fine- 
tuning condition leading to a set of physically interesting solutions. Since the geometries 
described in this paper are regular (i.e. curvature singularities are absent) gravity can be 
described in classical terms both in the bulk and on the vortex. We studied a thin string 
limit and identified the choice of parameters that may potentially lead to a solution to the 
gauge hierarchy problem. 

Various questions are still open. Our explicit solutions couple together scalar, tensor and 
gauge degrees of freedom. Therefore they represent an ideal framework where the localization 
of fields of various spin can be explicitly analyzed in a completely regular geometry. There 
are also open questions concerning higher windings. It has been shown that also for higher 
windings there are solutions which localize gravity. It would be natural to ask if these 
solutions are as stable as the ones obtained in the case of the lowest winding. Finally, since 
the geometries discussed in this paper are completely regular (in a technical sense) it is 
interesting to understand their stability against first order (quadratic) corrections to the 
Einstein-Hilbert action (which may naturally arise in a string theoretical context). 

We thank P. Tinyakov and K. Zuleta for discussions. This work was supported by the 
FNRS grants no. 21-55560.98, 7SUPJ62239 and by the Tomalla Foundation. 
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FIG. 21. A solution in the Bogomolnyi limit for c <C 1. 
APPENDIX A: THE BOGOMOLNYI LIMIT 

In this appendix we show analytically that for the case a = 2 (the Bogomolnyi limit 
the following limits hold: 



lim i/{c) = 2 

c— >0 

lim Mo (c) = 



c^O 



(A.l) 
(A.2) 



and 



lim Cn(c) 



(A.3) 



if the string tensions are tuned according to Eqs. ( [4.21[ )-( ^.22| ). The numerical solution for 



the Bogomolnyi limit is shown in Fig. 21 



In the Bogomolnyi limit and for c = the equations of motion for the gauge and scalar 
fields reduce to a coupled set of first order non-linear differential equations which can be 
written as 
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Cf' = fp, 

P' 



If we insert Eqs. (JKlj) into Eqs. (|]2TD-(|]2|) we obtain that 



Tp{x) = T0{x) = 0. 



The equation for M becomes then 



4m' + lOm^ = — /X. 



(A.4) 



(A.5) 



(A.6) 



The hmit c ^ imphes that /i ^ 0. From Eq. (|A.6|) the only (regular) solution respecting 
the boundary conditions at the origin is then M{x) = Mq = 1. Thus, this proves Eq. ( |A.2| ). 

Since Tp{x) = 0, from Eq. (^) we get 5^ = a/8 = 1/4 and from Eq. ([4.24| ) u = 2. 
Thus, this proves that, once the string tensions are tuned, Eq. ( |A.1| ) is valid. 

Now, the difference of Eqs. ( p.l3| ) and ( p. 141 ) yields. 



{M^M'L)' - {M^L'y = uM^L{fo - fe) ■ 



By using for P{x) Eq. ( p.l2| ) we have: 



PP'M^ 
aL 



(A.7) 



(A.8) 



Thus all terms in ( [A.7[ ) are total derivatives; integrating both sides from to x we obtain 

„ yPP' 



m 



a L2 



u P' 
l + -n — 
a L 



(A.9) 



where Mo is the value of M on the core of the string and n is the winding. In the case of 
fine-tuned solution, according to Eq. ( [4.24[ ), the last term in Eq. ( |A.9| ) vanishes. Since we 
just showed that in the limit a ^ 2 we have z/ ^ 2 and M — 1, from Eq. ( |A.9| ) we get 



LL' 



-PP' , 



which gives, after integration over x, 



L' = l- P' 



(A.IO) 



(A.ll) 



where we used the fact that -P(O) = 1. This proves ( |A.3|) , because P{x) — > for x — > oo. 
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